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Typhoid fever and salmonellosis, which are caused by Salmonella typhi and typhimurium, 
respectively, are responsible for significant morbidity and mortality in both developed 
and developing countries. We model typhoid fever using mice infected with Salmonella 
typhimurium, which results in a systemic disease, whereby the outcome of infection is 
variable in different inbred strains of mice. This model recapitulates several clinical aspects 
of the human disease and allows the study of the host response to Salmonella typhimurium 
infection in vivo. Previous work in our laboratory has identified three loci (Ity, Ity2, and 
Ity3) in the wild-derived MOLF/Ei mice influencing survival after infection with Salmonella 
typhimurium. Fine mapping of the Ity3 locus indicated that two sub-loci contribute col- 
lectively to the susceptibility of B6.MOLF-/fy//fy3 congenic mice to Salmonella infection. 
In the current paper, we provided further evidence supporting a role for Ncf2 (neutrophil 
cytosolic factor 2 a subunit of NADPH oxidase) as the gene underlying the Ity3. 7 sub-locus. 
Gene expression profiling indicated that the Ity3. 7 sub-locus defined a global gene expres- 
sion signature with networks articulated around Ncf2. Furthermore, based on differential 
expression and complementation analysis using Selp (selectin-P) knock-out mice, Selp was 
identified as a strong candidate gene for the \ty3.2 sub-locus. 
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INTRODUCTION 

Salmonella enterica, an intracellular Gram-negative bacterium, is 
the causative agent for a wide spectrum of clinical diseases with 
manifestations ranging from asymptomatic carriers, self-limiting 
gastroenteritis to fatal systemic infection (1, 2). There are over 
2500 Salmonella serovars, of which, some are host adapted such 
as serovar Typhi and Paratyphi in humans, while others, such as 
Typhimurium and Enteritidis, have a broad host range and are 
capable of infecting multiple organisms. In humans, Salmonella 
typhi causes a systemic disease, typhoid fever, which has a global 
health burden of 26.9 million cases and 200,000 deaths annually 
(3). In humans, Salmonella typhimurium is the causative agent of 
salmonellosis, a self-limiting gastroenteritis that results from the 
consumption of contaminated food or water. The emergence of 
multi-drug resistant strains of Salmonella in recent years high- 
lights the need for a more comprehensive understanding of the 
pathogenesis of Salmonella infection and for the identification of 
novel drug targets for vaccines and therapeutics (4). 

Salmonella is a natural pathogen of mice and infection with 
Salmonella typhimurium results in a typhoid-like systemic dis- 
ease. This murine experimental model has been used to iden- 
tify several genes and pathways involved in disease pathogenesis 
(5-9). As there is limited genetic variation within the classical 
inbred strains, the use of wild-derived strains of mice, such as 
MOLF/Ei contributes added genetic diversity and has allowed for 



the identification of novel genes that play an important role in 
innate immunity (10-13). Classical and wild-derived strains of 
mice exhibit a range of susceptibilities to Salmonella infection; 
for example, the C57BL/6J classical inbred strain are extremely 
susceptible to infection with Salmonella typhimurium due to a 
mutation in Slcllal ( solute carrier family 1 1 member 1 ) , while the 
129 sub-strains are highly resistant (14). The wild-derived mouse 
strain, MOLF/Ei is also susceptible to infection despite carrying 
functional copies of genes known to be important in Salmonella 
infection, such as Slcllal and Tlr4 (toll-like receptor 4) (6, 10). 

In order to identify the genetic determinants involved in 
the susceptibility of MOLF/Ei mice to Salmonella infection, 
we have previously used linkage analysis in an F2 panel of 
(C57BL/6 x MOLF/Ei) mice to identify two loci linked to host 
defense against Salmonella typhimurium, Ity2 (Immunity to 
Typhimurium locus 2) and Ity3 (10, 12). The MOLF/Ei allele at 
the Ity2 locus improves resistance to infection, whereas MOLF/Ei 
allele at the Ity 3 locus confers susceptibility (15). Validation and 
fine mapping of Ity3 locus were done using congenic B6.MOLF- 
Ity/Ity3 mice (12) and a panel of 12 sub-congenic mice (16). Using 
this approach, the Ity3 locus was refined to a 24 Mb interval and 
was shown to carry two sub-loci, Ity3.1 and Ity3.2 that together 
contribute to increased susceptibility to infection (16). The Ity3.1 
sub-locus controls NADPH oxidase activity and is characterized 
by decreased reactive oxygen species (ROS) production, reduced 
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inflammatory cytokine response, and increased bacterial burden. 
The Ity3.2 sub-locus is characterized by a hyper- responsive inflam- 
matory cytokine phenotype after exposure to Salmonella (16). 
Sequencing, expression, and functional data support the candi- 
dacy of Ncf2 (neutrophil cytosolic factor 2 a subunit of NADPH 
oxidase) as the gene underlying the Ity3.1 sub-locus (13). 

In the current study, we used global expression profiling to 
better understand the genetic networks that are being influenced 
by the Ity3 sub-loci and to identify potential candidate genes for 
the Ity3.2 sub-locus. We illustrate the impact of the Ity3.1 sub- 
locus on cell death and cytoskeletal reorganization, hematopoiesis 
as well as propose the candidacy of Selp (selectin P) as one of 
the candidate genes underlying Ity3.2 based on expression analy- 
sis, coding sequence polymorphism, and functional and allelic 
complementation studies. 

MATERIALS AND METHODS 
ETHICS STATEMENT 

All animals were maintained at the Animal Care Facility of McGill 
University according to the guidelines of the Canadian Council 
on Animal Care (CCAC). The animal protocol for this study was 
approved by the McGill University Animal Care Committee. 

ANIMALS 

Classical inbred strain C57BL/6J and wild-derived MOLF/Ei mice 
were used to generate congenic, B6.MOLF-ify and B6.MOLF- 
Ity/Ity3 and sub-congenic mice as described previously (12, 16). 
The susceptible Ity3 and resistant Ity mice, as well as the interme- 
diate B6MOL¥-Ity/Ity3.RecN and B6MOLV-Ity/Ity3.RecG mice 
were used for the microarray expression analysis, while the 
B6.MOLF-Jfy and B6.MOLF-.fty/7fy3, crossed with B6.129S7- 
Sei ,ptmiBay / j on j e red from the Jackson Lab (Bar Harbor, ME, USA), 
were used for the complementation assay. 

IN VIVO SALMONELLA INFECTION 

Mice aged 7-12 weeks were infected with Salmonella typhimurium 
strain Keller as described previously (12, 16). Briefly, mice were 
inoculated with 0.2 ml of physiological saline containing 10 3 
colony-forming units of bacteria through the caudal vein. The 
infectious dose was verified by serial dilutions on trypticase soy 
agar. Mice were either monitored for survival or euthanized at day 
3 or day 5 post-infection for organ collection. The animals were 
monitored two to three times daily and mice showing body con- 
dition scoring <2.0 were used for clinical endpoint (17). Survival 
analysis was conducted using a Kaplan-Meier survival test. 

MICROARRAY EXPRESSION ANALYSIS 

RNA was extracted from the spleens of mice, which were collected 
before infection and at day 3 post-infection. The RNA extrac- 
tion was carried out using TRIzol reagent (Invitrogen Canada, 
Inc., Burlington, ON, Canada). Three age-matched male mice 
were used per group. The concentration of RNA was determined 
using a NanoDrop spectrophotometer (Thermo-Fisher Scien- 
tific, Waltham, MA, USA). All hybridization and scanning of 
mice microarrays were carried out at the McGill University and 
Genome Quebec Innovation Centre, using the Illumina BeadArray 



technology (Illumnia Inc., San Diego, CA, USA). The expres- 
sion data were analyzed using FlexArray and normalized using 
a Lumi algorithm (Illumnia). Following the normalization, two 
approaches were used to generate a list of genes differentially 
expressed across the various strains. First, a Cyber f-test (Baladi 
and Long) (18) was used to generate a list of genes differen- 
tially expressed during infection, by comparing the fold change 
in expression of each gene between infected and uninfected sam- 
ples for each strain. These gene lists were further refined using the 
Benjamini Hochberg false discovery rate algorithm. Genes with an 
FDR p-value < 0. 1 and a fold change of > 2 were used to generate a 
final list of genes, which represented genes differentially regulated 
in each strain upon infection. 

A second approach was used to generate a list of genes that 
were differentially expressed at each time point, as compared to 
the control Ity strain. This was done by comparing the expression 
of genes in each strain to the control Ity strain at both day 0 and 
day 3 post-infection. The gene lists were further refined using the 
Benjamini Hochberg false discovery rate algorithm. Genes with an 
FDR p-value of <0.05 was used as a cut-off to characterize genes 
as significantly differentially regulated as compared to Ity. The 
gene lists generated using the two approaches were studied using 
a suite of online tools including DAVID (19), GeneGo (MetaCore, 
Thomson Reuters) and Gene Mania (20, 21). 

SEQUENCING OF Selp 

Sequencing was performed on PCR-amplified cDNA from Ity and 
Ity3 congenics obtained by reverse transcription of Trizol spleen- 
extracted RNA to determine genetic variation between C57BL/6J 
and MOLF/Ei alleles of the Selp candidate gene. Sanger sequenc- 
ing was completed at the McGill University and Genome Quebec 
Innovation Centre. 

ALLELIC COMPLEMENTATION ASSAY 

In order to study the effect of a MOLF/Ei Selp allele on suscepti- 
bility to infection, we carried out a complementation cross. Selp 
knock-out mice, B6.129S2-Se/p tmlHyn /J (Selp-'-) were ordered 
from the Jackson Laboratories (Bar Harbor, ME, USA). These 
mice were on C57BL/6J background with a mutant Slcllal allele. 
In order to correct for this, we crossed the Se/p~'~ mice to Ity 
as well as to the Ity3 mice. Mice were genotyped for the Se/p~'~ 
allele, and mice carrying the MOLF/Ei Slcllal allele along with the 
Selp~ ! ~ allele were inter-crossed to generate homozygous Selp~ ! ~ 
mice with a MOLF/Ei allele at the Slcl lal gene. Furthermore, these 
mice were crossed with Ity or Ity3 mice to generate mice, which are 
homozygous wild-type at Slcl 1 al, but carry a Selp knock-out allele 
complemented by either a C57BL/6J or MOLF/Ei allele (Figure SI 
in Supplementary Material) . 

BACTERIAL LOAD ENUMERATION 

For bacterial burden quantification in the spleen and the liver, mice 
were euthanized using CO2 and at the required day post-infection; 
both organs were removed aseptically, weighed and homogenized 
using a Polytron (Kinematica, Bohemia, NY, USA). The resulting 
homogenate was diluted in 0.9% saline and plated on tryptic soy 
agar to determine organ bacteria burden. 
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STATISTICAL ANALYSIS 

Statistical analysis was performed using Graph Pad Prism 6 
(GraphPad Software, San Diego, CA, USA). One-way ANOVA with 
Dunnet's multiple correction test was used to analyze the bacte- 
rial burden in the spleen. A corrected p-value <0.05 was used to 
establish significant differences. 

RESULTS 

Ity3 INFLUENCES THE EXPRESSION OF SPECIFIC GENES AND 
PATHWAYS DURING INFECTION 

The Ity3 locus is a complex QTL containing at least two sub-loci. 
We have previously studied the phenotypic impact of these two 
sub-loci (16), and in order to further characterize the genetic net- 
works that are affected by the MOLF/Ei allele at the Ity3 locus, we 
have used a microarray expression approach. Using this approach, 
we should be able to characterize and identify candidate genes 
within each sub-loci and identify pathways, which affect the host 
immune response to Salmonella infection. We studied the splenic 
global expression profile of infected and uninfected congenic resis- 
tant B6.MOLF-Iry (Ity), susceptible B6.MOLF-7ry/7fy3 (Ity3), and 
Ity3.RecG and Ity3.RecN strains. These two sub-congenic strains 
were selected because they carry either the MOLF/Ei alleles at 
Ity3.1 (Ity3.RecG) or at Ity3.2 sub-locus (Ity3.RecN), which result 
in an intermediate survival phenotype after infection with Salmo- 
nella typhimurium (16). Ity, Ity3, Ity3.RecN, and Ity3.RecG mice 
show significant differences in spleen and liver bacterial burden 
by day 5 post-infection (16). To evaluate primary defects in gene 
expression rather than differences due to high bacterial load, we 
selected day 3 as the time point to be studied because there was no 
significant difference in the spleen bacterial burden among the four 
strains of mice (16). Therefore, any changes in gene expression will 
serve as an indicator for genes and pathways that are important in 
regulating the susceptibility to Salmonella infection and not a con- 
sequence of a difference in bacterial load within strains. Our aim 
was to identify transcriptional signatures common to both Ity3 
and Ity3.RecG to define pathways controlled by the Ity3.1 locus 
and those similarly regulated by Ity3 and Ity3.RecN to identify 
transcriptional networks controlled by the Ity3.2 locus. We used 
two approaches to study the gene expression profiles (Figure 1). 

As an initial approach, we studied the changes in splenic gene 
expression in each strain over time (Day 3-Day 0). Overall, 241 
(Ity), 204 (Ity3), 218 (Ity3.RecG), and 201 (Ity3.RecN) genes were 
differentially expressed during infection as defined by a cut-off of 
a fold change >2 and a p-value <0.1 (Tables S1A-D in Supple- 
mentary Material; Figure 1). We identified 129 genes that were 
commonly and significantly regulated during infection in all four 
strains including a number of pro-inflammatory genes (Table S1E 
in Supplementary Material). The majority of these genes were 
chemokines, cytokines, and other immune related genes with a 
significant number of genes regulated by Type I and Type II 
IFN including several members of the Gbp family, Statl, Uspl8 
as well as others that are known to be involved in Salmonella 
infection (7, 15). These genes (Table S1E in Supplementary Mate- 
rial) define a transcriptional signature that is common during 
Salmonella infection and has been previously detected in different 
strains of mice during infection (15,22). Additional strain-specific 
genes regulated during infection were detected only in Ity (40 



genes), Ity3 (20 genes), Ity3.RecN (62 genes), and Ity3.RecG (54 
genes) (Tables S1F-I in Supplementary Material). A large num- 
ber of differentially expressed genes specific to Ity mice (Table 
S1F in Supplementary Material) were up-regulated in granu- 
locytes and/or macrophages including S100a8 and S100a9 that 
are of particular interest as they are involved with expression of 
inflammatory mediators, phagocytosis, oxidative burst as well as 
migration of neutrophils and monocytes to the site of infection 
(23). The gene Clec7a (dectin 1) was differentially regulated only 
in Ity. Recent work has linked dectin 1/Syk kinase signaling with 
autophagy-dependent maturation of phagosomes (24). 

In Ity3 mice, most genes that were differentially regulated were 
either expressed in macrophages or in megaerythrocyte precur- 
sors suggesting active extramedullary erythropoiesis in this strain 
(Table S1G in Supplementary Material). Interestingly, a large pro- 
portion (-40%) of genes specific to the strain Ity3.RecG are known 
to be down-regulated in B and T cells, as analyzed by BioGPS 
(25) (Table S1G in Supplementary Material). These data show the 
impact of Ity3 on the cellular composition of the spleen and/or 
changes in gene expression in specific cellular populations during 
infection. Very few genes were similarly regulated between Ity3 
and the sub-congenic strains Ity3.RecG (4 genes), Ity3.RecN (7 
genes), and Ity3.RecN and Ity3.RecG (4 genes) (Tables S1J-L in 
Supplementary Material). 

We then used GeneGO (Thomson Reuters, NY, USA) to clas- 
sify the genes differentially regulated in each strain (Tables S1A-D 
in Supplementary Material) into gene ontology (GO) molecu- 
lar pathways, GO processes, pathways, and process networks, in 
order to identify the pathways differentially regulated in each 
strain during infection (Figure 2). The strain Ity3.RecG appeared 
to have fewer genes involved in various pathways and processes 
related to chemokine and cytokine activity and immune response 
as demonstrated by the lower -log2 (p-values). These results are 
consistent with previous observations of reduced inflammatory 
responses following in vivo Salmonella infection in Ity3.RecG 
mice (16). 

Ity3.1 LOCUS INFLUENCES THE EXPRESSION OF GENES INVOLVED IN 
CELL-CYCLE REGULATION AND HEMAT0P0IESIS 

We have previously reported that Ncf2 is a strong candidate for the 
Ity 3.1 locus (13). To further characterize the downstream impact 
of Ity3.1 and Ncf2 on gene expression, we analyzed the data by 
evaluating variation in gene expression between Ity3, Ity3.RecG, 
or Ity3.RecN, and the control strain Ity at day 0 and day 3 post- 
infection (Figure IB). Interestingly, 12 genes located within the 
introgressed Ity3 region had lower expression levels in Ity3 and 
Ity3.RecG, compared to both Ity and Ity3.RecN. Figure 3A shows 
the box plot expression pattern of one of these genes and the 
list is provided in Table S2A in Supplementary Material. None 
of these genes were differentially regulated during infection in 
any of the sub-congenic strains. For some of these genes, the low 
levels of expression detected in mouse strains Ity3 and Ity3.RecG 
appear to be a consequence of poor hybridization of MOLF/Ei 
cDNA to microarray probes as a result of high genomic variability 
between the MOLF/Ei and C57BL/6J strains (Figure S2 in Supple- 
mentary Material). The wild-derived inbred MOLF/Ei had been 
separated from the classical inbred trains by over 1 million year 
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FIGURE 1 | Figure illustrating the two approaches used to study 
the changes in gene expression during infection of congenic and 
sub-congenic mice. (A) A representation of an initial approach used 
for studying the change in gene expression over time. The gene 
expression at day 3 post-infection was compared to day 0, and the 



of evolution, and as a result they have accumulated significant 
sequence variability, to the order of 1 SNP every 100 bp (26, 27). 

A second set of genes showed similar patterns of up-regulation 
(Figure 3B; Table S2B in Supplementary Material) or down- 
regulation (Figure 3C; Table S2C in Supplementary Material) dur- 
ing infection in all strains, although the constitutive and induced 
expressions were similar in Ity3 and Ity3.RecG but significantiy 
different from Ity and Ity3.RecN. This expression pattern high- 
lights the impact of the MOLF/Ei allele at the Ity3 locus, as all the 
genes, which an expression pattern similar to Figures 3B,C carry 
a MOLF/Ei allele at the Ity3.1 locus (16). 

Another set of genes showed higher expression during infec- 
tion only in strain Ity3 and Ity3.RecG (Figures 3D,E; Tables S2D,E 



final number of genes differentially regulated is stated on the right. 
(B)The second approach used, where gene expression in each strain 
was compared to the control Ity strain, at both day 0 and day 3. The 
final number of genes within each category is shown at the bottom of 
the panel. 



in Supplementary Material) as compared to Ity. In contrast 
to this grouping, the cluster of genes in Figure 3F exhibit 
lower expression levels in Ity3 and Ity3.RecG. Collectively, these 
groups of genes in Figures 3D-G (Tables S2D-G in Supple- 
mentary Material) exhibit a pattern of expression that is sim- 
ilar between Ity3 and Ity3.RecG, and the expression differ- 
ences could not be attributed solely to the MOLF/Ei allele at 
chromosome 1, therefore, we can conclude that the expres- 
sion differences are likely under the influence of the Ity3.1 
locus. 

Functional annotation of genes, which have an expression pat- 
tern similar between Ity3 and Ity3.RecG (Tables S2G-D in Sup- 
plementary Material) showed that a large percentage of the genes 
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FIGURE 2 | Pathway and process analyses of genes differentially 
regulated in the spleen of It y, Ity3, lty3.RecG, and lty3.RecN mice during 
infection. Gene ontology classifications generated by clustering the genes 
that are differentially regulated in each strain upon infection. (A) Gene 
ontology (GO) molecular pathways and (B) Gene ontology (GO) processes 
(C) pathways (D) Process networks that are enriched in the four mouse 



strains. The sub-congenic strain lty3.HecG shows a lower — log(p-value) for 
chemokine and cytokine receptor activity within the GO molecular pathways 
as well as lower — log(p-values) for response to bacteria in the GO processes. 
This is consistent with previous data suggesting that the lty3.RecG mice have 
a diminished inflammatory response following infection compared to parental 
strains (16). 



in the list plays a role in cell cycle, DNA binding, cytoskeletal 
reorganization, and hemopoietic and lymphoid organ develop- 
ment (Table S3 in Supplementary Material) (19, 28). Another 
major category of genes such as Ankl and Uros are involved in 
heme metabolic process. These data are consistent with the obser- 
vation that ROS control cell-cycle progression by influencing the 
presence and activity of cyclins and cyclin dependent kinases (29) 



and with a role for ROS in maturation and lifespan of erythroid 
cells (30,31). 

Ity3.1 AFFECTS A GROUP OF GENES THAT ARE CO-EXPRESSED WITH 
Ncf2 

In order to define the gene expression profile of the suscep- 
tible strains, we identified genes that had a similar pattern of 
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FIGURE 3 | Genes that are under the control of the It y3. 1 locus 

Sample box plots of the gene lists provided in Table S2 in Supplementary 
Material are shown. (A) Represents genes, which do not show any 
changes in expression during infection and show a similar expression 
pattern in Ity3 and lty3.FtecG. (B,C) Represents a sample box plot of 
gene expression, which show a similar regulation pattern in Ity3 and 
lty3.RecG. (D-G) Also illustrate genes withinTablesS2D,G in 



Supplementary Material, which show similar expression patterns in Ity3 
and lty3.RecG, but different from Ity and lty3.RecN. These genes are 
likely under the control of the Ity3. 1 locus. Csrpl (cysteine and 
glycine-rich protein 1), Chi3l1 (chitinase-like 1), (Plekhm3) Pleckstrin 
homology domain containing, family M, member 3), Sdsl (serine 
dehydratase-like), Cd59a (CD59a antigen), Lpxn (Ieupaxin), Asprvl 
(aspartic peptidase, retroviral-like 1). 



expression in the susceptible strain, Ity3 as well as the two sub- 
congenic strains Ity3.RecN and Ity3.RecG. Figures 4A,B highlights 
two genes, Tor3a and Fam20b as examples of the expression pat- 
tern of the list of genes provided in Table S2H in Supplementary 



Material, which have a similar expression pattern in Ity3, Ity3.RecN, 
and Ity3.RecG. Only 7 of the 47 genes were within the Ity3 interval, 
and almost all of them were within the genomic region common 
to Ity3.RecN and Ity3.RecG. This gene list was classified within 



Frontiers in Immunology | Microbial Immunology 



August 2014 | Volume 5 | Article 375 | 6 



Khan et al. 



Ity3 contribution during Salmonella infection 




www.frontiersin.org 



August 2014 | Volume 5 | Article 375 | 7 



Khan et a 



Ity3 contribution during Salmonella infection 



FIGURE 4 | Continued 

Genes that are differentially expressed in all susceptible and 
intermediate strains are under the influence of Ity3. 1 sub-locus. 

Box plots of the expression pattern of two genes (A) Tor3a (torsin family 
3, member A) and (B) Fam20b (family with sequence similarity 20, 
member B) are shown as examples to illustrate the expression pattern 
seen in the gene list provided inTable S2I in Supplementary Material. 
This list of genes show a similar expression in Ity3, lty3.RecG, and 
lty3.RecN and highlights the complex nature of the Ity3 locus as there 
are multiple genes in which expression is influenced by the combination 
of the two sub-loci. The genes that show a similar expression pattern in 



Ity3, lty3.RecG, and lty3.RecN were studied using GeneMania and the 
results are shown (CD), with the query genes being highlighted in 
black. Genes, which are known to be co-expressed, co-localized, have 
shared domains or predicted interactions with the list of genes in Table 
S2I in Supplementary Material are shown. The functional categories, 
which are enriched within this gene list, are inflammatory response, 
angiogenesis, and integrin mediated signaling pathways shown in 
yellow, blue, and red, respectively. Genes that were not differentially 
expressed but important in these pathways are shown in gray. 
(D) Co-expression of query genes, as well as other genes within these 
pathways, with Ncf2 is shown. 



functional categories (Figure 4C) such as inflammatory response 
and regulation of angiogenesis. A large proportion of the genes 
were either co-expressed, co-localized, or have shared domains or 
predicted interactions with Ncf2 (Figure 4D). We have previously 
shown that the MOLF/Ei allele at the Ity3. 1 locus contributed the 
strongest effect on susceptibility to Salmonella infection and was 
responsible for high bacterial burden and low ROS and cytokine 
production (16). The fact that a number of genes differentially reg- 
ulated in Ity3, Ity3.RecN, and Ity3.RecG strains, were co-expressed 
with Ncf2, supports the important contribution of the Ity3. 1 locus 
on the pathogenesis of infection in MOLF/Ei and its interaction 
with the other sub-locus Ity3.2 to enhance the impact of Ncf2 on 
ROS production. 

Selp IS A CANDIDATE GENE FOR Ity3.2 

We next studied genes showing a similar regulation pattern in 
Ity3 and Ity3.RecN to understand the pathways differentially regu- 
lated in Ity3.2 and identify potential candidate genes for the Ity3.2 
locus. There were only 14 genes that showed a similar expres- 
sion pattern in Ity3 and RecN (Figure 5). Of these 14 genes, 12 
are located on chromosome 1 and 6 genes (F5, Pbxl, Cacybp, 
Bc055342, Selp, and Vamp4) lie within the genomic region har- 
boring Ity3.2 (Table 1). Sequence variations have been reported 
between the MOLF/Ei and C57BL/6J in coagulation factor F5, the 
cDNA BC055324 and selectin P (Selp). The coagulation factor 
V is synthesized by the liver and is involved in the acceleration 
of prothrombin to thrombin conversion (32). Coagulation Fac- 
tor V deficiency leads to a bleeding disorder associated with mild 
to severe hemorrhagic symptoms (33). The cDNA BC055324 is 
poorly characterized and its function is not known. The Selp gene 
encodes for an adhesion molecule that mediates the recruitment 
of immune cells to the site of inflammation and is critical for the 
host immune response to infection making this gene an attractive 
candidate gene for Ity3.2. 

We further evaluated the candidacy of Selp as the gene under- 
lying Ity3.2 using sequence analysis and complementation assay 
in vivo. Selp encodes for a protein of 768 amino acids with a C-type 
lectin domain, an EGF domain and 8 complement control protein 
(CCP) modules [or as short consensus repeats (SCRs) functional 
domains] . We re-sequenced the coding region of Selp in C57BL/6J 
and MOLF/Ei mice and identified eight SNPs (Table 2), all of 
which are within the less homologous CCP domains, involved 
in protein recognition processes (34). The amino acid proline 
at position 205 is well conserved across 12 mammalian species 



and P205S is only observed in the DBA sub-strains, which share 
MOLF/Ei ancestry for this region of the genome (mouse phylogeny 
viewer) (35). In order to validate that the sequence variation in the 
MOLF/Ei Selp gene has an impact on susceptibility to infection 
with Salmonella typhimurium and to evaluate if Selp is indeed the 
gene underlying Ity3.2, we used an allelic complementation assay 
(see breeding scheme in Figure SI in Supplementary Material). Ity 
{Selp Mm ) and Ity3 (Selp MOLF/MOLF ) mice were crossed to Selp-'- 
knock-out mice and susceptibility to infection was assessed by 
survival analysis in Fl progeny with Selp MOLF/ ~ and Selp B6/ ~ 
genotypes. Selp M0LF/ ~ mice were significantly more susceptible 
to infection than Selp B6/ ~ mice and Ity controls (Figure 6A). We 
observed a lack of complementation in Selp 1 ^ 01 ^^ mice with 
a mean survival time equivalent to Selp~ j ~ animals (MST of 
8.1 ±0.18 and 8.3 ±0.38, respectively), adding further support 
for the candidacy of Selp as the gene underlying the Ity3.2 locus 
(Figure 6A). Although the SeZp~'~ mice showed a similar suscep- 
tibility compared to the Selp MOLFI ~ mice in terms of survival, their 

tissue bacterial burden was significantly lower when compared to 
SelpMOLF/-^ Jfy3> and RecG 

mice (Figure 6B) suggesting that the 
Ity3.2 locus does not contribute significantly to the bacterial bur- 
den and that the high bacterial burden observed in Selp MOLFI ~ 
mice is rather the effect of the Ity3.1 locus. 

DISCUSSION 

The current study was specifically designed to understand the 
pathways that are influenced by the Ity3 locus using sub-congenic 
strains that exhibit different degrees of susceptibility to Salmo- 
nella infection. The global gene expression profile in the spleen 
was studied early during infection prior to a significant bacterial 
difference in the target tissue. This approach allowed us to identify 
networks, which are of importance in the early phases of innate 
immunity yet not influenced by the extent of bacterial burden in 
the spleen. We reported a group of genes, the majority of which 
are regulated by type I and type II IFN. These genes, such as the 
Gbp, Oas, Ifitm family members, are differentially regulated in 
all strains of mice during infection, and define a core transcrip- 
tional signature common to several strains of mice infected with 
Salmonella (15). 

Additionally, we characterized a number of genes not located 
within the Ity3 region that were differentially expressed in the 
susceptible Ity3, Ity3.RecG, and Ity3.RecN strains as compared to 
the resistant Ity strain. We showed that these genes are also co- 
expressed with Ncf2, further supporting the hypothesis that there 
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FIGURE 5 | Continued 

Genes under the influence of Ity3.2 List of genes showing a similar 
expression pattern in Ity3 and lty3.RecN, and different from Ity and 
lty3.RecG. Ints7 (integrator complex subunit 7), F5 (coagulation factor V), 
Pbxl (pre B cell leukemia homeobox 1), Cacybp (calcyclin binding protein), 



Nenf (neuron derived neurotrophic factor), Angel2 [angel homolog 2 
(Drosophlla)], Cox6a2 (cytochrome c oxidase subunit Via polypeptide 2), 
Selp (selectin, platelet), Vamp4 (vesicle-associated membrane protein 4), 
Sccpdh [saccharopine dehydrogenase (putative)], Ctsw (cathepsin W), Fhl 
(fumarate hydratase 1). 
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FIGURE 6 | Allelic complementation between Ity3 congenic mouse and 
Selp deficient mice. In order to assess the impact of the MOLF/Ei Selp 
locus, Selp' 1 ' mice were crossed with Ity and Ity3 mice to generate mice 
carrying the knock-out allele complemented by the MOLF/Ei allele at Selp 
or C57BL/6J allele at the Ity3 locus. (A) Survival curves of the congenic Ity 
(n = 7), Ity3 (n = 7), knock-out Selp-'- (n = 21), and Se/p-' M0LF (n = 36), 
Selp- !m (n = 56) mice after infection with Salmonella typhimurium.Jhe 
Selp' 1 ' and Selp~ m0Lf show a similar curve after infection with Salmonella 
typhimurium, both of which are more susceptible than the control Ity 
congenic mice and Selp~ m , but more resistant than the mice carrying the 
entire Ity3 locus [Ity3). (B) Bacterial burden in the spleen of congenic Ity 
(n = 8) and Ity3 (n = 6), sub-congenic lty3.RecG (n = 8), and lty3.RecN 
(n = 8) and Se/p~'~ (n = 6) and compound heterozygous (n = 7 for Selp- /MoiF 
and n= 10 for Selp~ IBB mice after infection). 



Table 1 | List of known SNPs within differentially expressed genes in 
Ity3 and lty3.RecN 

Target Variation Chromosome Location (bp) 



DOWN-REGULATED 



Ints7 No exonic variation 




1 


191575734 


F5 rs6271495 




1 


164151838 


Pbxl No exonic variation 




1 


168153527 


Cacybp No exonic variation 




1 


160202367 


D730003i15rik No exonic variation 




1 


191224474 


Nenf No exonic variation 




1 


191306789 


Bc055324 rs30651611 




1 


163945993 


Angel2 No exonic variation 




1 


190925112 


Cox6a2 No exonic variation 




7 


128205436 


UP-REGULATED 








Selp rs30667849 




1 


164115264 


Fh1 rs 13465421 




1 


175600374 


Vamp4 No exonic variation 




1 


162570515 


Sccpdh No exonic variation 




1 


179668210 


Ctsw No exonic variation 




9 


5465240 


Table based on data from the Mouse Genome Informatics (MGI) website. 


Table 2 | Exonic variation in the MOLF/Ei allele of the Selp gene. 


Base pair change Amino acid change 


Domain 


Selp 730G >T 


V202F 




Sushi/CCP/SCR 


739A > G 


N203D 




Sushi/CCP/SCR 


745C >T 


P205S 




Sushi/CCP/SCR 


620C>T 


H207Y 




Sushi/CCP/SCR 


841G>A 


G239S 




Sushi/CCP/SCR 


1135G>A 


V337I 




Sushi/CCP/SCR 


1775A>C 


N550T 




Sushi/CCP/SCR 


1831A>G 


I569V 




Sushi/CCP/SCR 



is an influence of the Ity3. 1 sub-locus on other segments of the 
genome. These results highlight the importance of the region of 
distal chromosome 1 carrying Ity3, a region enriched in QTLs. 
Over 80 QTLs are listed at the mouse genome database (36) and 
a number of cis and trans eQTLs have been characterized (37). 
Among them, several QTLs are involved with complex inflamma- 
tory reaction, such as graft vs. host disease (38) lupus (39), modifier 
of LPS-response (40), and susceptibility to tuberculosis (41). 

We have reported previously that Ity3 influences ROS produc- 
tion during infection (13) and this effect was mapped recently to 
a small sub-region named Ity3.1, which harbors the gene Ncf2, a 
subunit of the NADPH complex. (16). ROS produced by NADPH 
has been shown to affect a number of pathways, which are impor- 
tant in innate immunity including bacterial killing within the 



phagolysosome where ROS interact with other ions such as chlo- 
ride to form toxic agents (HOCL) or can convert into hydroxyl 
radicals that are toxic for bacteria (42). ROS production has also 
been shown to influence immune cell recruitment, activation, 
and survival (43), and lead to translational activation of NF- 
kB (44). NADPH oxidase derived ROS is also a key regulator 
of autophagy and autophagy regulation during pathogen inva- 
sion (45). In addition, NADPH oxidase activation contributes to 
the recognition and removal of apoptotic neutrophils (efferocyto- 
sis) by macrophages (46-48). Therefore, an imbalance in NADPH 
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produced ROS can lead to increased inflammation, which can be 
deleterious to the host. 

Recent studies have shown that TNF, as well as other cytokines 
through NF-kB signaling induces transient increase in ROS level 
in endothelial cells, which results in cell surface expression of Icam 
and Selp (49-52). This effect has further been studied in vivo, where 
mice lacking the p47 phox (Ncfl) subunit of the NADPH complex 
have reduced expression of VCAM-1, ICAM-1, SELR and SELE in 
the vascular cell walls (53, 54). In the current study, we identified 
Selp as a candidate gene for the locus Ity3.2 and used an allelic 
complementation assay to provide genetic evidences that Selp is 
indeed a strong candidate for the Ity3.2 locus. In MOLF/Ei mice 
susceptibility to infection as explained by the Ity3 locus could 
be attributed to the individual effect of Ity3.1 and Ity3.2 sub- 
loci but also to the cooperation between these two sub-loci as 
explained by the potential impact of low activity of NADPH oxi- 
dase (16) on Selp function by reducing its expression (current 
paper). 

We also illustrated that mice carrying a MOLF/Ei allele at the 
Ity3. 1 locus have higher expression of a number of genes playing a 
role in cell cycle, DNA binding, and cytoskeletal pathways. There 
is a growing body of evidence discussing the link between ROS, 
cell-cycle progress and arrest. As discussed by Martindale (55), 
ROS can have multiple effects on the cell cycle, depending on the 
amount and type of ROS. They suggest that low doses of ROS may 
cause proliferation while high doses of ROS can lead to apoptosis 
and cell death. In our model, it is possible that reduced levels of 
ROS in susceptible animals could lead to cell growth arrest, hence 
providing a more favorable niche for bacteria to replicate. 

Another pathway influenced by ROS production, is up- 
regulated by Ity3. 1 and regroups genes involved in heme biosyn- 
thesis. Increased expression of genes within the heme biosynthesis 
pathway could result in increased free heme within the cells, which 
can act as a potent cytotoxic pro-oxidant (56). Free heme has also 
recently been shown to trigger necroptosis in macrophages (57). 
Therefore, it is possible that the increased expression of the heme 
biosynthesis pathway observed in susceptible mice is a mechanism 
that compensates for the low ROS levels. 

In conclusion, our study highlights the complex nature of 
multi-loci interaction in the wild-derived MOLF/Ei response to 
Salmonella infection. We highlighted the role of low ROS and 
cytokine production in reduced survival of mice carrying the 
Ity3.1 locus, and the importance of the Ity3.2 locus, which syn- 
ergistically led to increased susceptibly of the lty3 mice. We 
have also shown that several pathways identified in strains Ity3, 
Ity3.RecN, and Ity3.RecG, are influenced by Ncf2. Furthermore, the 
Ity3.1 sub-locus has additional effects, which have not previously 
been characterized, in expression of genes involved in cell-cycle 
arrest and hematopoiesis. Lastly, we propose a hypothesis that 
the combined effects of low ROS production by the MOLF/Ei 
Ity3.1 locus together with the impact of Selp MOLF/Ei allele at 
Ity3.2 influences the host survival after infection with Salmonella 
typhimurium. 
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